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Abstract. The third dredge- up phenomenon in asymp- 
totic giant branch (AGB) stars is analyzed through evo- 
lutionary model calculations of a 3 M Q , solar metallicity 
star. The Schwarzschild criterion is used to test the stabil- 
ity of a given layer against convection, and the calculations 
are performed either with or without extra-mixing below 
the convective envelope. Based on these calculations, sev- 
eral questions are addressed regarding the occurrence of 
the third dredge-up in AGB star models, the laws govern- 
ing that phenomenon, and some of its implications on the 
structural and chemical evolution of those stars. 

The use of the Schwarzschild criterion without extra- 
mixing of any sort is shown to lead to unphysical afterpulse 
models which prevent the occurrence of third dredge-up. 
Model calculations of a 3M Q star using no extra-mixing 
confirm the failure to obtain dredge-up in those condi- 
tions. That conclusion is found to be independent of the 
mixing length parameter, stellar mass, or numerical accu- 
racy of the models. 

Model calculations performed on selected afterpulses of 
the 3 M star, but with extra-mixing (using a decreasing 
bubble velocity field in the radiative layers and a diffusion 
algorithm for the mixing of the chemical elements), lead to 
efficient dredge-ups at a rate of 10~ 5 — 10~ 4 M Q /yr. Test 
calculations using different extra-mixing extents and effi- 
ciencies reveal that the dredge-up predictions are rather 
insensitive to those extra-mixing parameters. This impor- 
tant conclusion is understood by analyzing the physics 
involved in the dredge-up process. It is shown that the 
dredge-up rate is determined by the thermal relaxation 
time-scale of the envelope as C-rich matter is added from 
the core into the envelope. The dredge-up predictions are, 
however, expected to depend on the convection prescrip- 
tion in the envelope. 

Linear relations both between the dredge-up rate and 
the core mass M c and between the dredge-up efficiency A 
and M c are predicted by the model calculations. Those lin- 
ear relations are expected to still hold when the feedback 
of the dredge-ups on the AGB evolution is taken into ac- 
count. They predict the dredge- up efficiency to level off at 
unity during the AGB evolution, at which point the core 
mass remains constant from one pulse to the next. The 



core mass is concomitantly predicted to evolve towards 
an asymptotic value. The existence of such an asymptotic 
core mass naturally provides an upper limit to the mass 
of the white dwarf remnant, and helps to constrain the 
initial-final mass of white dwarfs. 

Synthetic calculations taking into account the dredge- 
up laws obtained from the full AGB model calculations 
predict a continuous increase of the stellar luminosity L 
with time, contrary to the predicted behavior of M c and 
A. This results from an adopted dependence of L on both 
M c and the radius R c of the H-depleted core of the form 
L oc M^/Rc- As a result of this increase of L with time, 
the initial-final mass relation can further be constrained 
if mass loss is taken into account. If, for example, a su- 
perwind is assumed to eject all the remaining envelope 
of the 3Mq star at L = 15000 L Q , then the mass of the 
white dwarf remnant is predicted to be 0.66 M Q , instead 
of 0.73 Mq predicted by models without dredge-up. 

Finally, the synthetic calculations predict the forma- 
tion of a 3 Mq carbon star after about 20 pulses experi- 
encing dredge-up. Taking into account the fact that the 
luminosity decreases by a factor of two during about 20% 
of the interpulse phases, such a 3 M Q carbon star could be 
observed at luminosities as low as 7500 L Q . 

Key words: stars: AGB - stars: carbon - stars: white 
dwarf - stars: evolution - stars: interior - stars: abundances 



1. Introduction 

The third dredge-up (3DUP) phenomenon is the process 
by which ashes from helium burning nucleosynthesis are 
transported from the interior to the surface of low- and 
intermediate-mass stars (M <, 7M Q ) during the asymp- 
totic giant branch (AGB) phase (see, e.g., Mowlavi 1998a 
for a general introduction to the structural and chemical 
evolution of AGB stars). It contrasts with the first and sec- 
ond dredge-ups which bring to the surface only ashes from 
hydrogen burning. The signature of those 3DUP events is 
observable in the specific chemical composition character- 
izing many AGB stars as compared to those of red giant 
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stars (see, for example, Mowlavi 1998b). In particular, the 
high C/O ratios observed in MS, S, SC (C/O <, 1 in num- 
ber) and C (C/0>1) stars are explained as resulting from 
an increase in the carbon abundance due to the 3DUP 
(solar C/O~0.42). Similarly, the overabundances (as com- 
pared to those in red giants) of fluorine and s-process el- 
ements (Jorissen et al. 1992) attest to the operation of 
dredge-up in AGB stars. 

Stellar models show that thermal instabilities (called 
pulses) develop periodically in the He-burning shell of 
AGB stars (Schwarzschild and Harm 1965), leading to 
structural readjustments which provide an adequate stage 
for the operation of the 3DUP (Iben 1975, Sugimoto and 
Nomoto 1975). Unfortunately, model predictions available 
in the literature do not agree on the characteristics of the 
3DUP, such as the pulse number at which it begins to 
operate along the AGB or the extent of envelope pene- 
tration in the C-rich layers. The problem is of particular 
relevance when confronting chemical abundances observed 
at the surface of AGB stars to model predictions. For ex- 
ample, the famous 'carbon star mystery' put forward by 
Iben (1981) about the luminosity function of C stars in 
the Magellanic Clouds revealed, among other things, the 
difficulty of models to explain the low luminosity C stars 
(see Mowlavi 1998c for the present state on that ques- 
tion). The few - but isolated - successes of AGB models 
in obtaining low-luminosity carbon star models in the late 
80 's suggested that the solution to the carbon star mys- 
tery resides in the use of new radiative opacities and ade- 
quate mixing length parameters (Sackmann & Boothroyd 
1991). Yet, many recent AGB calculations do not confirm 
those expectations (Vassiliadis & Wood 1993, Wagenhu- 
ber & Weiss 1994, Blocker 1995, Forestini & Charbonnel 
1997), while some others do succeed in obtaining dredge- 
up (Frost & Lattanzio 1996, Straniero et al. 1997, Herwig 
et al. 1997, this work). 

The disagreement between AGB model predictions on 
the issue of dredge-up prevents us to provide reliable and 
consistent chemical yields of low- and intermediate-mass 
stars. Why do some models obtain dredge-up while oth- 
ers fail to do so? Herwig et al. (1997) succeed in obtain- 
ing dredge-up from the 7 th pulse on in a 3M© Pop. I 
star model by using an overshooting algorithm below the 
convective envelope. Some sort of extra-mixing beyond 
the convection borders defined by the Schwarzschild crite- 
rion was already used by Boothroyd & Sackmann (1988a) 
and Lattanzio (1986), and its inclusion is found to im- 
prove the dredge-up efficiencies predicted by the models 
of Frost & Lattanzio (1996). These last authors further 
note the critical dependency of dredge-up on the numer- 
ical treatment of convection in massive AGB models. In 
contrast, Straniero et al. (1997) claim to obtain dredge- 
up consistently without invoking any extra-mixing. A sim- 
ilar claim is also reported by Lattanzio (1989). Accord- 
ing to Straniero et al. (1997), the solution resides rather 
in the time and mass resolution of the models. Clearly, 



several questions must be clarified. How influential is the 
numerical accuracy of AGB models in obtaining dredge- 
up? What is the role of extra-mixing below the envelope? 
Is extra-mixing necessary for obtaining 3DUP in current 
AGB models using the Schwarzschild criterion? If yes, how 
sensitive are the dredge-up characteristics to the extra- 
mixing parameters, about which little is known today? 
And, finally, what are the characteristics governing the 
3DUPs? 

The aim of this paper is to shed some light on these 
questions. They are, for a great part, related to the deli- 
cate question of convection and the determination of their 
boundaries in AGB models. The use of a local theory 
of convection, such as that of the mixing length the- 
ory (MLT), is certainly one of the greatest shortcomings 
still affecting stellar model calculations in general, and 
AGB models in particular. New non-local formulations 
are being developed (e.g. Canuto & Dubovikov 1998) and 
may, hopefully, be included in future calculations. For the 
present time, however, evolutionary AGB models are still 
performed using the MLT prescription - mainly due to 
computer time requirements -, and we also adopt it in 
this study. During the third dredge-up, however, the use 
of such a local prescription leads to the development of an 
(unphysical) discontinuity in the abundance profiles when 
the H-rich envelope penetrates into the H-depleted layers. 
The determination of the lower boundary of the convec- 
tive envelope then becomes problematic, and requires a 
careful analysis of its stability against mixing across the 
discontinuity (i.e. extra-mixing). That problematic is set 
forth in Sect. ||. Several definitions involving the use of 
the Schwarzschild criterion and of extra-mixing are also 
presented in that section. Sections [| to ^ then analyze 
the 3DUP predictions from AGB model calculations of a 
3 Mq star of solar metallicity, computed both with and 
without extra-mixing. Section ^ describes the main char- 
acteristics of the code used to compute the models. A first 
set of calculations performed without any extra-mixing is 
analyzed in Sect. || It is shown that dredge-up does not 
occur in those conditions. In contrast, Sect. ^| shows that 
models calculated with extra-mixing lead naturally to the 
occurrence of the 3DUP phenomenon. The dredge-up laws 
derived from those calculations are presented in Sect. [| 
Finally, some implications of those laws on stellar evolu- 
tion along the AGB are analyzed in Sect. 0. Conclusions 
are drawn in Sect. |[ 



2. On the convective boundaries in AGB stars 

The criterion used in most stellar evolution calculations - 
and in all AGB models published to date - to delimit con- 
vective regions is that of Schwarzschild. According to that 
criterion, a layer is convective if V r ad > V a d, and radiative 



N. Mowlavi: On the third dredge-up phenomenon in asymptotic giant branch stars 



3 



if V ra d < V Q d- V a d is the adiabatic temperature gradient, 
and V ra d the radiative temperature gradient given by 



in model stars where both of the above cases are present 



simultaneously. This is discussed in Sect. 2.3. 
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with T, L, k, and m being, respectively, the tempera- 
ture, luminosity, Rosseland mean opacity, and mass con- 
tained in the sphere interior to the radius r (other symbols 
having their usual meanings). The layer where V, a d = 
V a d defines the convection boundary, and is called the 
Schwarzschild layer. 

Let us suppose for now that mixing of material inside 
a star can only occur through convection, and let us imag- 
ine the star as being composed of a finite number of shells 
located at increasing mass coordinates to. Each shell is 
characterized by a temperature T(m,t), density p(m,t), 
luminosity L(m, t) and chemical composition X(to, t). All 
these quantities vary with time t. The Schwarzschild cri- 
terion tests locally the stability of each shell against con- 
vection by comparing V ra d{rn, t) derived from Eq. |l| with 
V a d at the given shell independently of the values of those 
quantities in adjacent shells. Mixing of matter between 
two adjacent shells is allowed to operate only if both of 
them arc found to be convective. Alternatively, the chem- 
ical composition of a radiative shell is not altered with 
time, according to the Schwarzschild criterion, unless and 
until its radiative gradient exceeds its adiabatic gradient. 
This could result from a change with time of T, p and/or 
L (we neglect here changes due to nuclear burning). Ap- 
plied in this strict definition, the local Schwarzschild crite- 
rion is called hereafter the strict Schwarzschild criterion. 
We note that the shell-representation of a star used in 
the above discussion corresponds to the numerical zon- 
ing performed during the computation of stellar models 
(by the finite difference method). Thus, stellar model cal- 
culations using the Schwarzschild criterion without any 
extra-mixing essentially use the strict Schwarzschild crite- 
rion. For this reason, the 'strict Schwarzschild criterion' is 
equivalently called, hereafter, the Schwarzschild criterion 
without extra-mixing. 

Convection, however, is essentially non-local. Convec- 
tive bubbles reaching the convection boundary do pene- 
trate, on a whatever small distance, into the radiative lay- 
ers. The chemical composition of a radiative layer adjacent 
to a convective zone may thus, in reality, be modified by 
this convective penetration. For it to occur in models us- 
ing the Schwarzschild criterion, however, an extra-mixing 
procedure must be specified in the code (whether purely 
numerically or according to a physical prescription). In 
most cases, this extra-mixing does not affect the stability 
of the radiative layers adjacent to the convective zone. A 



short discussion of these cases is presented in Sect. 2.1 



In some cases, however, this extra-mixing does modify the 
stabil ity of those radiative layers. Those are discussed in 
Sect. 2.2. They characterize, in particular, the dredge-up 
phase in AGB stars. A third situation may be encountered 



2.1. Case a: Stable Schwarzschild boundary 

The most common case is encountered when the convec- 
tion boundary lies in a region of smooth V ra( z profile, 
such as in chemically homogeneous regions as illustrated 
in Fig.gfaJ. This case characterizes, in particular, the bot- 
tom layers of the convective envelope in AGB stars during 
most of their pulse and interpulse phases. 

It is easily seen that an extra-mixing of material from 
the convective envelope into the radiative layers does not 
alter the stability of those layers against convection. The 
Schwarzschild layer is thus said to be stable in the sense 
that an extra-mixing of matter from the convective zone 
into the underlying radiative layers does not, to first order, 
alter the location of the Schwarzschild layer. 

Of course, the Schwarzschild layer does not necessarily 
locate correctly the boundary of the convective zone, since 
a penetration of the convective bubbles into the radiative 
layers (i.e. overshooting) may occur, the extent of which is 
still a question of debate. But such an overshooting does 
not, to first approximation, have any crucial consequence 
on the location of the Schwarzschild layer. 

2.2. Case b: Unstable Schwarzschild boundary 

During a 3DUP, the H-rich convective envelope penetrates 
into the H-depleted radiative layers, and gives rise to a 
discontinuous chemical profile [dotted line in Fig. [TJfb )] . As 
a result, a concomitant discontinuity develops in the V ra£ j 
profile [solid line in Fig. ^(b); see also Iben 1976, Paczynski 
1977, Frost & Lattanzio 1996]. A similar situation occurs 
during the core helium burning phase of most stars when 
the He-depleted convective core grows into the overlying 
He-rich layers. 

In models using the Schwarzschild criterion without 
extra-mixing, the border of the convective zone is auto- 
matically located at the position of the hydrogen abun- 
dance discontinuity. As a result, no alteration of any chem- 
ical abundance is expected to occur in the radiative layers. 
The convection boundary remains located at this position 
as long as the value of W ra d on the radiative side of the 
discontinuity is below W a d- Evolutionary AGB model cal- 
culations performed without extra-mixing confirm that no 
penetration of the convective envelope occurs beyond that 
point of discontinuity (see Sect. ||). It should be stressed 
that the Schwarzschild layer is undefined in those models, 
a situation which is clearly unphysical. 

In reality, some matter does travel across the V ra d 
discontinuity. This results from the finite velocity of the 
bubbles at the convection boundary, since V ra d — V a£ £ is 
positive. Models should thus include some sort of extra- 
mixing. The consequences of such a penetration can read- 
ily be evaluated: mixing of matter from the envelope into 
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Fig. 1. Hydrogen mass fraction (dotted lines), adiabatic 
temperature gradient (dashed lines) and radiative tem- 
perature gradient (solid lines) profiles at the bottom of 
the convective envelope in three different AGB stars as 
a function of their mass coordinate: (a) during the shell 
H-burning phase of a 3M star; (b) during the afterpulsc 
phase of a 3M Q star; (c) during the afterpulse phase of a 
6M Q star. 



the radiative regions enriches the latter region with hydro- 
gen and renders them unstable against convection. This in 
turn leads to further penetration of the envelope. The sit- 
uation is thus unstable, in the sense that an extra-mixing 
of matter from the convective zone into the underlying 
layers does alter the location of the convective boundary. 

The challenge facing AGB modelers is to provide a re- 
liable description of this envelope penetration with simple 
(i.e. computer feasible) convection prescription. The anal- 
ysis presented in Sect. || actually reveals that the penetra- 
tion of the envelope into the H-depleted layers is rather 
insensitive to the extra-mixing parameters such as its ex- 
tent or mixing efficiency. That main conclusion to be ob- 
tained in Sect. ^| is most welcomed since not much - not 
to say nothing - is known on the extra-mixing character- 



istics in AGB stars. Section || also shows that the use of 
a diffusive algorithm for the mixing of chemical elements 
enables to avoid the development of a V r ad discontinuity 
in the models, and to properly define the Schwarzschild 
layer. 

2.3. Case c: Existence of a minimum in the V r ad profile 

In massive AGB stars, the V ra( j profile can display a mini- 
mum close to the lower boundary of the envelope, as illus- 
trated in Figjfjfc). A situation can occur where two convec- 
tive zones develop, separated by a radiative region around 
that minimum. This case combines the features of both 
previous cases: the outer convective zone has a smooth 
V ra d profile as in case a, while the inner one presents a 
discontinuity as in case b. 

In models using the Schwarzschild criterion with no 
extra-mixing, the inner convective zone does not penetrate 
into the H-depleted regions due to the V ra d discontinuity 



(see Sect. 2.2). The outer zone, on the other hand, is free to 
deepen until it reaches the inner convective zone. The two 
convective zones then merge together, and the situation 
resumes to case b. 

When extra-mixing is allowed to operate, a specific 
scheme is expected to develop. Because of the small 
amount of mass involved in the inner convective zone 
(~ 10 _4 M Q ), any extra-mixing into the He-rich layers 
reduces the hydrogen content of that convective zone, and 
renders its outermost layers stable. Those outermost layers 
are then progressively excluded from the convective zone 
as they reach convective neutrality, and a semi-convective 
zone is expected to develop]]. Eventually, the main convec- 
tive envelope penetrates into that semi-convective zone, 
and the case resumes to case b. 

2.4- Conclusions 

It is clear from the discussion in this section that mod- 
els using the Schwarzschild criterion with no extra-mixing 
are inadequate to describe the 3DUP phenomenon. There 
is, however, a great deal of confusion in the literature on 
the third dredge-up predictions, and in particular on the 
role of extra-mixing. For the purpose of clarifying some of 
those issues, AGB models have been computed both with 
and without extra-mixing. Particular attention has been 
devoted to ensure that the code comply with the defmi- 



A semi-convective zone is defined here to be one throughout 
which partial mixing leads to convective neutrality V rat j = V ' a d 
(see, e.g., Mowlavi & Forestini 1994). This definition is not to 
be confused with other ones used in the literature, such as that 
whereby a semi-convective region is defined by V a£ j < V ra d < 
V ' Led, V 'Led being the Ledoux adiabatic gradient (which takes 
into account the gradient of the mean molecular weight). Fur- 
thermore, the semi-convective zone discussed here, and located 
at the H-He discontinuity, is different from that described by 
Iben (1981), which is located at the He-C discontinuity. 
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tions given at the beginning of this section on the location 
of convection boundaries. The description and results of 
those calculations are presented in Sects. || to ||. 



3. The stellar evolution code 

The stellar evolution code used to compute the models is 
described in detail in Mowlavi (1995). Among the modifi- 
cations brought to it since its conception in 1995, let me 
mention the use of interior radiative opacities taken from 
Iglesias & Rogers (1996) and the use of low-temperature 
ones from Alexander & Ferguson (1994). The convection 
energy transport is described as usually by the MLT, with 
the mixing length parameter a m it taken equal to 1.5 in 
units of pressure scale height. The time-step and mesh al- 
location procedures rel evan t to our study are presented 
in more details in Sect. 3A, while the mixin g schemes of 
chemical elements are described in Sect. 3.2. 



3.1. Time- step and mesh distribution 

The time step between each model is determined, in most 
cases, by the relative variation of the dependent variables 
(radius, temperature, pressure and luminosity) from one 
model to the next such that it does not exceed 10% at 
any layer in the star. On the AGB phase, the time step 
is further limited to a maximum of A^ rater . p /1000, where 
Atinterp is the time between two successive pulses. The 
time-step during the 19th interpulse phase, for example, 
is of the order of 25 y, and can decrease down to few hours 
during a pulse. 

The mesh allocation is based on the constitutive dif- 
ferential equations, and ensures a relative accuracy in the 
finite difference method^ of better than 5 x 1CP 5 (see 
Mowlavi 1995). Special care is taken around boundaries 
of convective regions and in regions with steep chemical 
abundance profiles in order to avoid spurious numerical 
chemical diffusion. Our models are characterized by about 
3000 meshes during the interpulse phases, and by up to 
8000 meshes during pulses. 

3.2. Mixing 

Instantaneous mixing is assumed in all convective zones 
in models with no extra-mixing. When a discontinuity in 
the Vrad profile is detected at the convection border, care 
is taken not to mix the material of the mesh comprising 
that discontinuity with the convective zone. This prevents 
a spurious numerical propagation of the convective zone 
into the radiative layers in unstable cases such as case b 



Our criterion for the mesh allocation contrasts with that 
used in almost all existing stellar evolution codes where the 
mesh distribution is determined according to the relative vari- 
ation of the dependent variables (except that of Wagenhuber 
& Weiss 1994 who use a method similar to ours). 



described in Sect. 2. 2, and guarantees a correct application 
of the Schwarzschild criterion described in Sect. |[ 

In models with extra-mixing, a specific overshooting 
prescription is used beyond the Schwarzschild layer. The 
bubble velocity field v(r) in the radiative field is as- 
sumed to decay exponentially with the distance to the 
Schwarzschild layer, in the form 



v(r) = v sch x x(r) x exp 
with 



MU) x 



(r sc h - r) 



(j"sch Tov^ 



fv 



V 8c h 



(2) 



(3) 



where subscripts sch and ov refer to quantities evaluated 
at the Schwarzschild boundary r sc h and at the edge of 
the overshooting region r ov , respectively. The function 
x{r) = \{r— r ov )/(0.5 r+0.5 r sc h — r ov )] 0A ensures that v(r) 
vanishes at the outer edge of the overshooting region. It is 
close to unity in all the overshooting region except close 
to r v where it vanishes. Mixing in the convective + over- 
shooting regions is performed through an algorithm cou- 
pling diffusive mixing and nucleosynthesis (see Mowlavi 
1995) with a diffusion coefficient equal to D conv (r) — 
\xv(r)xl mix (r), where l mix {r) = a m i t xH p (r) (H p being 
the pressure scale height). The bubble speed at r sc h, v sc h, 
is determined from the bubble velocity profile in the con- 
vective zone as given by the MLT, while that at r OVl v OVl 
is determined by an 'efficiency' parameter f v according to 
Eq. |^. The extent of the overshooting region (r sc h — f v) is 
determined by a second parameter, a ov , defined such that 
r ov locates the layer where P ov = Psch X exp(a 01J ), P ov 
and P sc h being the pressure at r ov and r sc h, respectively. 

In those models with extra-mixing, the mesh distribu- 
tion is increased in the overshooting regions by about 100 
meshes in order to ensure a good description of the abun- 
dance profiles. The time-step is further reduced to about 
5 x 10 -2 y during a 3DUP. A dredge-up episode is then 
covered by about 2000-3000 models. 

Let us remark that a velocity field of the type of Eq. || 
is supported by recent hydrodynamical simulations of stel- 
lar convection in solar-type stars and white dwarfs (Frey- 
tag et al. 1996). However, we do not claim that such an 
overshooting pattern also applies to the envelopes of AGB 
stars, in which the physical conditions are quite different 
from those characterizing solar-type stars or white dwarfs. 
That remains to be confirmed by future hydrodynamical 
studies. We rather use such a prescription merely as a way 
to simulate an extra-mixing below the envelope and to 
analyze its effects on the 3DUP predictions. The extra- 
mixing could very well be associated with other 'non- 
standard' mixing processes such as, for example, diffusion 
induced by rotation or shear mixing. Actually, the results 
presented in Sect. ^| show that the dredge-up predictions 
are quite insensitive to the extra-mixing parameters. 
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4. Models without extra-mixing 

4--1- Failure to obtain dredge-up 

A set of standarcQ 3 M models of solar metallicity is com- 
puted with a m it = 1.5 from the pre-main sequence phase 
up to the 32nd pulse on the TP-AGB phase. The maxi- 
mum depth reached by the convective envelope after each 
pulse is shown in Fig. |[ The depth is expressed in terms of 
the mass separating the lower convective envelope bound- 
ary to either the H-He discontinuity (lower solid line) or 
the He-C one (upper solid lines). It is seen that the en- 
velope does not penetrate into the He-rich layers (and a 
fortiori into the C-rich layers) during the first thirty two 
pulses. 

The non-occurrence of the third dredge-up in those 
models is easily understood as resulting from the dis- 
continuity in the H abundance when the convective en- 
velope reaches the H-He discontinuity (case b described 
in Sect. 2.2). In those regions, the opacity is mainly de- 
termined by the Thompson opacity Kth = 0.2(1 + X), 
where X is the hydrogen mass fraction. It presents a jump 
across the composition discontinuity, with the H-rich lay- 
ers being more opaque to radiation than the He-rich layers 
[Fig. ^(b)]. This discontinuity in k translates into a similar 
discontinuity in V ra d [Fig. 0(aj], which prevents further 
penetration of the envelope into the H-depleted regions 
when no extra-mixing is allowed in the models. 

It is instructive to recall that the structural readjust- 
ments characterizing the afterpulse phases result from the 
evacuation of a gravothermal energy (defined by e gra v+e = 
+ if , where P, p and e are the pressure, den- 



d(l/p) 
dt 



sity and internal energy, respectively) which develops at 
the bottom of the former pulse (Paczyhski 1977). This 
translates into a positive gravothermal energy wave prop- 
agating outwards and reaching the H-rich layers in about 
250 yr during the 19th afterpulse in our 3M Q models 
(see Fig. ^). The luminosity at the bottom of the enve- 
lope, L env> b, increases concomitantly with time until the 
gravothermal energy has been evacuated [Fig. As 
a result, V ra <2 increases too, and the convective envelope 
penetrates inwards until it reaches the H-He discontinuity. 
A further deepening of the envelope could be possible if 
the increase in L enVt b overcomes the discontinuity in V ra£ j. 
In the standard 3 Mq models, however, this does not hap- 
pen during the first thirty two pulses calculated so far. 

Effect of mixing length parameter The mixing length pa- 
rameter is known to be very influential on the depth of 
convective envelopes (Wood 1981). In order to analyze its 
role in obtaining 3DUP in AGB models without extra- 
mixing, the standard 3M Q star is recalculated from the 
first to the 23rd pulse with «mlt = 2.2 (i.e. similar to the 
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3 The word 'standard' refers throughout this paper to models 
calculated without any extra-mixing beyond the Schwarzschild 
layer. 



Fig. 2. Depth reached, as a function of the core mass, 
by the envelope of 3M Q , 2=0.02 models (where Z is the 
metallicity) during the first 20 to 30 thermal pulses (iden- 
tified by the filled circles and squares on each curve) . The 
two lower curves, labeled M r ^ env — M r ^ e , indicate the 
mass, in solar mass units, separating the lower boundary 
of the envelope to the location of the H-He discontinuity. 
The two upper curves give a similar information, but to 
the location of the He-C discontinuity. Solid and dashed 
lines refer to models calculated with olmlt = 1.5 and 2.2, 
respectively. 



value used by Straniero et al. 1997). The results are shown 
in dotted lines in Fig. 0. Indeed the convective envelope 
reaches deeper layers with olmlt = 2.2 than with 1.5. 
But they are not more successful in penetrating the He- 
rich layers. This contrasts with the results of Straniero et 
al. who claim to obtain dredge-up in their models without 
any extra-mixing. 

Effect of stellar mass Finally, a last set of standard mod- 
els is calculated for a 6M Q star of solar metallicity up 
to its 22nd pulse. It is indeed known that 3DUP should 
also be favored in more massive AGB stars (Wood 1981). 
The results of our model calculations are shown in Fig. || 
Again, no dredge-up is found in those models calculated 
without any extra-mixing, for reasons similar to those put 
forward for the 3 M Q standard models. 

Effect of numerical accuracy The analysis presented 
in this section shows that models using the local 
Schwarzschild criterion without extra-mixing do not lead 
to the occurrence of third dredge-up because of the V ra[ j 
discontinuity at the core's edge. That conclusion should 
thus not be sensitive to the numerical accuracy of the 
models. Indeed, test calculations performed on the 19th 
afterpulse without extra-mixing, but with increased or de- 
creased accuracies on both the time-step and mesh reso- 
lution lead to results identical to those presented in this 
section. 
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Fig. 3. (a) Radiative gradient profile at several times, as 
labeled (in yr) on the curves, after the 19th pulse of the 
standard 3 M case. The time origin is arbitrary but iden- 
tical to that of Fig. ^. The dashed dotted line is the adia- 
batic gradient taken at the time corresponding to the solid 
line; (b) Same as (a), but for the opacity profile. The ex- 
tents of the H-, He- and C-rich regions are shown at the 
bottom of (b), separated by vertical lines 



£ gra v + e ( 1 8 x erg/g *s) (a) 
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Fig. 4. (a) Gravothermal energy profiles at 5 different 
times, as labeled (in yr) on the curves, after the 19th pulse 
of the standard 3M©, solar metallicity star. The time ori- 
gin is taken at maximum pulse extension. The He-rich 
region is identified by the thick vertical line separating 
the H- and C-rich regions; (b) Same as (a), but for the 
luminosity profiles. The thick dots on the curves indicate 
the location of the lower convective envelope boundary 
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Fig. 5. Same as Fig. g but for 6 M , Z=0.02 models with 
CtMLT = 1.5 



4-2. Useful relations 

Some relations characterizing the evolution of the stan- 
dard 3M models are presented in this section. These 
will be useful in Sect. ||. 

The first relation describes the evolution of the surface 
luminosity L with core mass growth. The values of those 
quantities before the onset of each pulse of our standard 
3M Q models are displayed as filled circles in Fig. ^(a). It 
is well known that AGB models not experiencing dredge- 
up reach an asymptotic regime characterized by a linear 
relation between L and M c . This is the famous M c — L 
relation, which writes from our 3 M© models [dotted line 
in Fig. §(aj] 



L = 63340 (M c - 0.493), 



(4) 



L and M c being given in solar units. A correction has to 
be applied to this relation in order to account for the lower 
luminosities of the first pulses. This correction is found to 
be well reproduced by an exponential function, and Eq. ^ 
becomes 



L = 63340 (M c - 0.493) — 1100 x e -( 100A/c+58 ). 



(5) 



This relation is shown in solid line in Fig. ^(a), and is seen 
to fit very well model predictions. 

The second useful relation expresses the interpulse du- 
ration Ati n ter &s a function of core mass growth. From 
our standard 3M Q model calculations, shown in filled cir- 
cles in Fig. ^b), we find an asymptotic M c — log Ati nter p 
relation given by [dotted line in Fig. ^(b)] 



log At 



interp 



4.2 (1.8 -M c ), 



(6) 



Ati n terp being expressed in yr. Again, a correction has to 
be applied during the first pulses. The resulting relation 
writes [solid line in Fig. Wb)] 



log At 



interp 



4.2 (1.8- M c ) -0.1 x e-( 90M - 53 - 2 ). (7) 



Finally, the location R en v,b of the lower envelope 
boundary at its first minimum after the pulse is displayed 
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Fig. 6. (a) Stellar luminosities at maximum pulse exten- 
sions, (b) interpulse pulse duration between two maximum 
pulse extensions and (c) location of the lower envelope 
boundary at its first minimum after each pulse in the stan- 
dard 3 M Q model calculations (filled circles in each panel). 
The dashed lines show the linear relations fitting the data 
in each panel (Eqs. ||, || and ||). The solid lines in panels 
(a) and (b) show the M c — L (Eq. |J) an d M c — Ati nte rp 
(Eq. [?]) relations, respectively, when the deviation from 
linearity of the first pulses is taken into account. 



in filled circles in Fig. ||(c) as a function of M c . A linear 
relation fits the data in the asymptotic regime [dotted line 
in Fig. M(c)], which writes in solar units 



R 



env,b 



1.263 x (0.885 - M c ). 



(8) 



5. Dredge-up in models with extra-mixing 

5.1. Model results 

The structural evolution of the 19th afterpulse of the stan- 
dard 3M Q star is recalculated in a set of models using 
the overshooting prescription described in Sect. 3.2, with 
(ct ov ,f v ) — (0.60, 10 -4 ). The structural evolution of the 
intershell layers of these models is shown in long dashed 
line in Fig. ^(a). 

Clearly, the use of extra-mixing leads to the pene- 
tration of the convective envelope into the H-depleted 
regions, as expected from the discussion in Sect. 2.2. 
The role of overshooting in the 3DUP was already ad- 
dressed by Iben (1976) and Paczyhski (1977). The lat- 
ter author, in particular, showed the essential role of 
extra-mixing in obtaining 3DUP. The envelope is seen in 
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Fig. 7. Structural evolution of several quantities during 
and after the 19th pulse of the 3M star in the standard 
case (solid line) and in models with envelope overshoot- 
ing (dotted line: a ov = 0.20, f v = 10 -6 ; long dashed line: 
a ov = 0.60, /„ = 10 -4 ; short dashed line: a ov = 1.20, f v = 
10 -3 ). (a) Mass location of the convective boundaries 
and of the layer of maximum nuclear energy production 
(dashed-dotted lines) in the H and He burning shells, as 
labeled in the graph. Hatched regions correspond to con- 
vective zones. The extent of the overshooting regions is 
less than 10 _4 M© and cannot be distinguished in the fig- 
ure; (b) surface 12 C mass fraction; (c) stellar luminosity; 
(d) stellar radius; (e) location of the lower boundary of 
the envelope. The time origin is taken at maximum pulse 
extension. 

our models to deepen at a maximum dredge-up rate of 
4.2 x 10~ 5 M o /yr. The dredge-up lasts for about 120 yr, 
which corresponds to the time necessary for the gravother- 
mal energy to evacuate from the bottom of the formal 
pulse to the stellar surface (see Fig. ||). 3.6 x 10~ 3 M Q of 
C-rich material is mixed into the envelope, which leads to 
a 16% increase in the surface carbon mass fraction. 

It is instructive to stress that the extent of the over- 
shooting region below the envelope comprises less than 
10 _4 M Q . This is 40 times less than the total amount of 
material dredged-up. Overshooting alone can thus not ex- 
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Table 1. Several quantities characterizing the afterpulse phases of the 12th, 15th, 19th, 25th and 32nd pulses of the 
3 M star computed with the overshooting parameters a ov — 0.60 and f v = 10~ 4 : the core mass M c and the luminosity 
L at maximum pulse extension, the luminosity dip L m i n during the afterpulse phase [see Fig. Rfc)], the maximum 
dredge-up rate of the envelope rd U p, the core mass increase AM c nuc due to nuclear burning during the preceding 
interpulse phase, the mass of C-rich material dredged-up Md U p, the efficiency parameter A, the carbon mass fraction 
X p ( 12 C) in the intershell layers left over by the pulse, and the increase in the envelope 12 C mass fraction AX S ( 12 C) 
resulting from each dredge-up. 
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would have stopped from further penetration for reasons 
described in Sect. |]. 
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Fig. 8. (a) Same as Fig. Q(a), but for models with over- 
shooting (a ov = 0.60, /„ = 10~ 4 ). (b), (c) and (d) same 
as (a), but for the opacity, the gravothermal energy and 
the radius, respectively. 



plain the amount of material dredged-up. Indeed, the oc- 
currence of dredge-up essentially results from the unstable 
character of the convection border in the standard models 



(see Sect. 2.2). Extra-mixing has the role of triggering the 
penetration of the envelope into the H-rich layers. Had we, 
at any time during the dredge-up, suppressed the extra- 
mixing procedure in the model calculations, the envelope 



The use of a diffusive extra-mixing procedure re- 
sults in models displaying smooth \7 ra d profiles as 
shown in Fig. [1(a). This enables to locate properly the 
Schwarzschild layer. It should be noted, however, that the 
Schwarzschild layer is still unstable in the sense that more 
mixing of hydrogen from the envelope into the deeper 
layers renders them unstable against convection, leading 
thereby to further penetration of the envelope. 

The question then arises of the sensitivity of the 
dredge-up characteristics to the extra-mixing parameters. 
In order to answer that question, two other sets of mod- 
els are calculated, one with increased extra-mixing extent 
and efficiency of (a ov ,f v ) — (1.20, 10 -3 ), and one with 
lower extra- mixing extent and efficiency of (a ov ,f v ) = 
(0.20, 10~ 6 ). The predictions from those sets are displayed 
in Fig. in dotted and short-dashed lines, respectively. It 
is seen that the dredge-up characteristics are rather insen- 
sitive to the extra-mixing parameters. In particular, the 
dredge-up rate is independent of the extra-mixing param- 
eters, and equals 4.2 x 10~ 5 M Q /yr in all three sets of cal- 
culations. The only small difference between the different 
sets is the earlier time, in models with higher a ov , at which 
the envelope reaches the H-depleted regions. 

The important conclusion that the dredge-up charac- 
teristics are rather insensitive to the extra-mixing param- 
eters can be understood from basic considerations on the 
dredge-up process, developed in Sect. [| It allows, in par- 
ticular, to analyze the 3DUP phenomenon without wor- 
rying much about the exact values of the overshooting 
parameters. In order to perform such an analysis as a 
function of stellar parameters, additional calculations with 
extra-mixing are performed on the 12th, 15th, 25th and 
32nd afterpulses of the standard 3 M© star. The resulting 
structural evolution of the intershell layers are shown in 



Fig. I? 
Table 



and some of their characteristics summarized in 
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Fig. 9. Same as Fig. ^(a), but for, respectively from top 
to bottom, the 15th, 19th, 25th and 32nd pulses of the 
3M Q models with a ov = 0.60 and f v = 10~ 4 . The size of 
the mass ordinate is similar in all four figures. 

Dredge-up rate The dredge-up rate ra V p is seen from Ta- 
ble |l| to be an increasing function of the pulse number. 
A linear relation between r^up and M c is found from the 
model calculations, as shown in Fig. |l(](a). It translates 
into the following equation for our 3 M© star [dashed line 
in Fig. 

(9) 

in M Q and M yr^ 1 , re- 



r dup = 6 xlQ- 4 (M c - 0.595), 

M c and r^ U p being expressed 
spectively. 



Dredge-up efficiency The depth reached in mass coordi- 
nates by the envelope after each pulse depends on both the 
time over which it penetrates and the dredge-up rate. The 
former is mainly determined by the thermal time-scale for 
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Fig. 10. (a) Dredge-up rate and (b) dredge-up efficiency 
as a function of core mass in the models with a ov — 0.60 
and f v = 10~ 4 (filled dots). The dotted lines represent 
a linear fit to the data (Eqs. ^| and [l(] for Vd U p and A, 
respectively). The crosses and plus signs in panel (b) are 
the data predicted by Straniero et al. (1997) and Herwig 
et al. (1997), respectively. 



the evacuation of the gravothermal energy accumulated at 
the bottom of the pulse (Fig. while the latter is deter- 
mined by the time-scale of the thermal readjustment of the 
envelope. The resulting mass of C-rich material dredged- 
up in the envelope, Md up , is reported in Table [l| 

The dredge-up efficiency is traditionally expressed by 
a number A, which is the ratio of M du p to the core mass in- 
crease AM Ct „ uc due to hydrogen burning during two suc- 
cessive pulses. A value of A = 1, for example, implies a 
constant core mass from one pulse to the next. From the 
model calculations, A is found to also increase linearly with 
M c [filled circles in Fig. [l(](b )} . The linear relation writes 
[dotted line in Fig. ^(b)] 



A = 7 x (M c - 0.571), 

M c being expressed in solar mass. 



(10) 



Evolution of surface characteristics during the dredge-up 
The evolution of the surface luminosity and radius dur- 
ing the 3DUP is shown in Fig. J^(c,d) (solid lines), and 
compared to that predicted by models without dredge-up 
(dotted lines). It is seen that dredge-up leads to surface 
luminosities lower than those obtained in otherwise sim- 
ilar models but without dredge-up. That the luminosity 
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should be lower in models with dredge-up is easily under- 
stood from the above discussion, since the penetration of 
the envelope absorbs energy in order to lift the dredged- 
up material into the envelope (see also Paczyhski 1977). 
Dredge-up also leads to lower surface radii. 

5. 2. Influence of convection prescription 

While the dredge-up characteristics are found to be rather 
insensitive to the extra-mixing parameters, a dependency 
of those characteristics on the convection prescription in 
the envelope should be discussed. Higher values of a m i t 
are known to favor deeper convective envelopes (Wood 
1981). The envelope would thus reach the H-depleted lay- 
ers sooner and lead, for an otherwise similar dredge-up 
rate, to more C-rich material being dredged-up. Moreover, 
the dredge-up rate itself should also be affected by the con- 
vection prescription. The analysis presented in Sect. |6.1| 
stresses on the role of the thermal response of the envelope 
in the dredge- up process. The thermal time-scale of the en- 
velope depends on its thermal structure, which is imposed 
by the surface conditions and the a m u parameter. We thus 
expect the dredge-up rate to be also dependent on a m u- 
A test calculation on the 19th afterpulse of the 3 M com- 
puted with a m it = 2.2 leads to rd up = 6.9 x 1CP 5 M /yr 
with M C =O.65M . This has to be compared with a pre- 
dicted dredge-up rate of 5.7 x 10~ 5 M /yr at that core 
mass in models with a m it = 1.5 (from Eq. ||). Models with 
Oimit = 2.2 thus predict higher dredge-up rates by about 
16% than those in models with a m it = 1.5. 

5.3. Comparison with other published works 

The calculations presented in Sects. |] and [5] lead to a sim- 
ple conclusion about AGB models using the Schwarzschild 
criterion: models computed without any extra-mixing do 
not lead to dredge-up while those using an extra-mixing 
procedure lead to efficient dredge-up. 

How does this compare to model calculations published 
in the literature? The answer to that question requires 
the knowledge of the precise numerical procedures used 
by each AGB modeler in their code. Such information is 
often not available in the literature. For this reason - and 
because it is anyway outside the scope of this paper to 
do so - it is impossible to compare all published models. 
Some of them are however discussed below. 

Many AGB calculations of the last five years do not 
report the occurrence of the 3DUP (Vassiliadis & Wood 
1993, Wagenhuber & Weiss 1994, Blocker 1995, Forestini 
& Charbonnel 1997). This is consistent with their working 
hypothesis of the Schwarzschild criterion without extra- 
mixing. 

Several codes use an extra-mixing procedure based on 
a purely numerical technique. Boothroyd and Sackmann 
(1988a), for example, mix the radiative layer adjacent to 
the convective zone before checking for its stability. If the 



Vrad/Vad ratio of that (formally stable) layer is found 
to be greater than one after mixing, then they declare it 
convective. As this mixing procedure is performed once 
per model calculation, it allows a numerical propagation 
of the convective envelope during the 3DUP. This might 
be the reason for their success in obtaining low-mass car- 
bon stars (Boothroyd & Sackmann 1988b). A similar tech- 
nique was already used by Paczyhski (1977). This last au- 
thor obtained dredge-up in 8 M AGB models with such 
a numerical 'overshoot' procedure, while no dredge- up re- 
sulted in models not including that overshoot. The use of 
extra-mixing procedures based on a purely numerical tech- 
nique suffers from at least two shortcomings. First, they 
assume instantaneous mixing of the layer(s) added to the 
envelope. Besides being unphysical, it may lead to conver- 
gence difficulties. Second, they would probably depend on 
the numerical accuracies of the models such as the time- 
step and mesh re-zoning. As a result, such models may or 
may not lead to efficient dredge-up. 

Herwig et al. (1997) use an overshooting algorithm 
which assumes an exponentially decreasing bubble veloc- 
ity field in the extra-mixing region. Their algorithm is very 
similar to mine, and enables a direct comparison of their 
predictions with mine. These are reported in Fig. |lO|(b) 
(+ signs). Interestingly, they find higher dredge- up effi- 
ciencies than I do at a given stellar luminosity. We note 
that they use a m it = 1-7 in their calculations. This could 
explain part of the differences. We further note that their 
core masses are lower by about 0.01 M than mine at a 
given pulse number, which could imply a dependence of 
the dredge-up efficiencies on the initial conditions at the 
onset of the AGB. 

Finally, the recent work of Straniero et al. (1997) 
should be mentioned. These authors obtain dredge-up in 
models of 1.5 and 3M stars without using any extra- 
mixing procedure in their calculation. This seems in con- 
tradiction with my conclusions. Private discussions with 
some of the authors in Straniero et al. (1997) seemed 
to confirm that they did not use any extra-mixing in 
their model calculations. They explain their results by the 
higher temporal and spatial resolution of their models (see 
their paper for more details). This conclusion, too, is not 
supported by my calculations reported in Sect. || While 
I do not understand their results, I cannot push the anal- 
ysis much further. In any case, the dredge-up efficiencies 
found in their models are lower than those predicted from 
my calculations including extra-mixing. Their predictions 
are reported in Fig. ^Ojfbj (x signs). 

6. Third dredge-up laws 

6.1. The dredge-up process 

The important conclusion that the dredge-up rate is insen- 
sitive to the extra-mixing parameters can be understood 
from the following considerations. The AGB star is ba- 
sically formed by an e~-degenerate core of 0.10-0.15 R 
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surrounded by an extended H-rich envelope. As dredge-up 
proceeds, the outermost H-depleted layers of the core be- 
come part of the envelope and expand [which translates 
into a negative gravothermal energy production in those 
layers, see Fig. ^(c)]. This results in an increase of the po- 
tential energy of those layers engulfed in the convective 
envelope, which must be supplied by the luminosity pro- 
vided by the inner layers. The maximum dredge-up rate 
is then determined by the thermal relaxation time-scale of 
those layers expanding into the envelope during the dredge- 
up. 

Let us estimate that dredge-up rate from simplified ar- 
guments. The potential energy E pot of a mass Am located 
at the H-He discontinuity is given by 



E, 



pot 



-G 



Am Mr 



Br 



(11) 



where M c is the mass of the H-depleted core and R c its 
radius (i.e. R c — R en v,b during dredge-up). The energy 
necessary to lift the dredged-up material into the envelope 
is provided by the luminosity at the core edge, L c . The 
thermal time-scale tdup for the dredge-up is then 



td 



E, 



pot 



lip 



= 3.12 x 10' 



Rr L r 



A? 



(12) 



All quantities in this section are expressed in solar units 
and years. The dredge-up rate can then be estimated by 



T dup 



Am 

tdup 

= 3.2 x 10" 



Rc L c 

Mr 



(13) 



Typical values for the 3M star are M c ~ 0.66 M Q 
[Fig. |(aj], R c ~ 0.12 R [Fig. |(dj] and L c ~ 
10 4 L Q [Fig. ||(b)]. Equation |l3| then gives r^up — 
5.8 x 10~ 5 M Q /yr. This estimate of the dredge-up rate, 
while giving only an order of magnitude, convincingly sup- 
ports the rate of 4.2 x lO _5 M /yr obtained in the full 
evolutionary models (Sect. ||). It is also consistent with 
the results obtained by Iben (1976), who estimated the 
dredge-up rate to be about lO~ 5 M /yr in a 7M star 
with M c = 0.95 — 0.96 M Q , and with those of Paczyhski 
(1977), who found 5.2 x lO~ 5 M /yr in a 8M star with 
M c = 0.8 M . 

The question of why the dredge-up rate is indepen- 
dent of the extra-mixing parameters can further be under- 
stood by comparing the time-scale of the thermal read- 
justment of the envelope with the typical time-scale of 
convective bubbles to cross the H-He transition zone. The 
velocity of the convective bubbles as they approach the 
H-He discontinuity is of the order of 10Ro/yr, which 
translates in terms of mass to about lO _1 M /yr (from 
Fig. ||(d)). This rate is much higher than the dredge- up 
rate of ~ 10~ 5 M Q /yr established above. The deepening of 



the envelope into the H-depleted layers is thus primarily 
fixed by the thermal relaxation time-scale of the envelope, 
rather than by the speed of the convective bubbles pene- 
trating into the H-depleted layers. 

6.2. Dredge-up characteristics 

Equation [l^ reveals a formal dependence of the dredge-up 
rate on R c , M c and L c . 

Those three stellar parameters are not independent of 
each other. Refsdal & Weigert (1970) and Kippenhahn 
(1981), for example, show from homology considerations 
that most red giant star properties are functions of the 
H-depleted core mass and radius. As the burning shell 
advances, the structure of the intershell layers evolves, 
to first approximation, like homologous transformations. 
Let us consider two times t and t' characterized by core 
masses, radii and luminosities of M c , R c and L, respec- 
tively, at time t, and M' c , R' c and V , respectively, at 
time t'. Then the homology transformations applied to 
a ~ 3M AGB star lead to (Kippenhahn 1981, Herwig et 
al. 1998) 

2 / r,/ \ -1 

(14) 



V_ M^y R, 
L K \M C ) \R, 

The dredge-up rate r' dup 
uated from that at time t 
to 

r 'dup M' 



at time t' can then easily be eval- 
Tduv Equations [l^ and |l4]lead 

(15) 



Tdup 

The homology transformations thus suggest a linear 
relation between the dredge-up rate and the core mass. 

In standard AGB calculations (i.e. without dredge- up), 
R c is a linear function of M c (Eq. || for our 3M star). 
Equation |TJ then recovers the classical linear M c -L rela- 
tion. In that linear function of either M c , 
R c or L. The dredge-up predictions presented in Sect. || 
essentially obey these rules (Table [[]), since the dredge- up 
calculations have been performed on selected afterpulses 
of the standard 3 M models and do not include the feed- 
back of the dredge-ups on the AGB evolution. 

When the feedback of dredge-up on the structural evo- 
lution is taken into account, the core radius is no more a 
linear relation of the core mass, and the linear M c -L does 
not hold any more (Herwig et al. 1998). It must be re- 
placed by a M c -R c -L relation as suggested by Eq. [l4| 
Equation however, reveals that the dredge-up rate 
keeps a linear dependency on M c even in models expe- 
riencing dredge-up. 

The implications of the linear M c -rd U p relation on the 



core mass evolution are analyzed in Sect. 7.1. 



7. Dredge-up and the formation of white dwarfs 

The dredge-up laws established in Sects. || and || allow 
an approximate study of the structural evolution of AGB 
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stars without actually performing full AGB model calcu- 
lations. This is the so called 'synthetic' AGB calculation 
technique. 

Such a study is outside the scope of the present paper. 
However, a preliminary analysis of the synthetic evolution 
of a 3 M Q star with dredge-up is presented in Sects. ^ and 
7.3 using simplified assumptions. But before doing that, 
the implications of the linear dependence of the dredge-up 
rate on the core mass is first discussed in Sect. 

7.1. Core mass evolution and the formation of white 
dwarfs 

The analysis presented in Sect, ^supports a linear relation 
between rdup and M c . I assume that such a linear relation 
remains valid between A and M c even when the feedback 
from the dredge-ups on the structural AGB evolution is 
taken into account. 

Such a linear relation between A and M c has an impor- 
tant consequence on the core mass reached at the end of 
the AGB evolution. Equation [lC] predicts that A reaches 
unity at M* = 0.71 M©. It is easy to see that the feedback 
between A and M c leads to the evolution of M c towards an 
asymptotic value M c = M* . Indeed, if M c < M* , then A < 1 
and the core mass increases from one pulse to the next. 
If M c > M* , then A > 1 and the core mass decreases from 
one pulse to the next. At M c = M*, the amount of ma- 
terial dredged-up from the core during the 3DUP equals 
that added to the core by H-burning during the interpulse 
phase, and M c remains constant from one pulse to the 
next. We thus expect A to level off at unity, and M c to 
reach M*. This puts an upper limit to the mass of white 
dwarf remnants (which is 0.71 M Q for our 3M star), and 
helps to constrain the predicted initial-final mass relation 
for white dwarfs. 

7.2. Structural evolution 

According to Eq. [| dredge-up begins to operate in our 
3M Q star around the 10th pulse, at M c = 0.599 M . This 
is where we begin our synthetic evolution. We assume that 
the star has reached its asymptotic regime, which is a good 
enough assumption (see Fig. |^) for our purposes. 

The evolution of M c is governed, on the one hand, by 
its increase AM C „ UC due to hydrogen burning during the 
interpulse. For a solar metallicity star, AM Ci „ uc is given by 
(all quantities in this section are expressed in solar units 
and in years) 



AM Ci „„ c = /x, x 1.37 x 10 11 L H x At interp 
— fh x 1-37 x 10~ n L x At interpi 



(16) 



where the luminosity Lh of the H-burning shell is ap- 
proximated by the surface luminosity (at the time of the 
maximum extension of the next pulse). The factor f^ is 
introduced in order to account for the fact that the mean 
surface luminosity during the interpulse is actually lower 
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Fig. 11. Predictions from synthetic calculations for a 
3 M Q AGB star with (filled circles) and without (filled tri- 
angles) dredge-up. (a) Dredge-up efficiency, (b) core mass, 
and (c) surface luminosity as a function of time. The first 
displayed pulse (at time t = 0) corresponds to the 10th 
pulse of the standard 3 M Q models. The horizontal dashed 
line in (b) shows the asymptotic core mass predicted for 
models experiencing dredge-up. 

than that at the the next pulse. From the standard model 
calculations, we find fh = 0.72. After each pulse, on the 
other hand, M c decreases as a result of the 3DUP. The am- 
plitude of this decrease is given by Md up = A x AM c nuc , 
where A is estimated from Eq. The resulting net in- 
crease in M c from one pulse to the next is then 



AM C = (1 - A) X AM C ,, 



(17) 



In model calculations without dredge-up, the change in 
luminosity AL from one pulse to the next in the asymp- 
totic regime is given by Eq. 0, while the change in the core 
radius AR C is given by Eq. the presence of dredge-up, 
those relations must be modified. For AL, a dependence 
on (AM C ) 2 / AR C is suggested from Eq. |l4], and we adopt 



AL = -80000 



(AM C ) 2 
AR C 



(18) 



For Ai? c , we make the assumption that the time evolution 
of R c remains unaffected by the dredge-ups. Equations || 
and |^ then determine R c as a function of time t. 

Finally, the interpulse duration is assumed not to be 
altered by the dredge-up episodes, and is given by Eq. [?]. 

The predictions of our synthetic calculations are dis- 
played in Figs, [y] (time evolution) and [l^ (M c -L evo- 
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Fig. 12. Luminosity as a function of core mass predicted 
by the synthetic calculations for a 3 M© star when dredge- 
up is (filled circles) or is not (filled triangles) taken into 
account. 



lution). As expected from the discussion in Sect. 7.1, A 
increases towards unity and the core mass towards the 
asymptotic value of M* = 0.71 M . The luminosity, how- 
ever, increases continuously. This results from the fact 
that L depends not only on M c but also on R c according 
to Eq. [l4|. Had we used a linear M c — L relation, then L 
would have evolved towards an asymptotic value (of about 
14000 L Q for our 3M Q star) in a similar way as does A. 
The non-linearity of L with M c is made clear in Fig. |l2|. 

The assumptions underlying those synthetic calcula- 
tions certainly are too simplistic. For example, the feed- 



backs of dredge-up on A, At 



i nter p 



or AR r have been ne- 



glected, and relations [l4| and |l5| need to be confirmed by 
model calculations following the dredge-ups all along the 
AGB evolution. Yet, the main conclusions, such as the 
existence of a limiting M* towards which the core mass 
evolves asymptotically, should qualitatively be correct. 



Effects of mass loss A last word about the effects of mass- 
loss. The amplitude of mass loss depends mainly on the 
surface radius R and effective temperature T e f / (and thus 
on L). Those are expected to be modified by the dredge- 
ups (Fig. ^). However, the time evolution of L [Fig. Q(cj] 
reveals that they are not as much affected as is M c . As a 
result, mass loss is expected to be higher in models with 
dredge-up than in those without, at a given M c . In other 
words, mass loss will eject the AGB's envelope at a much 
lower value of M c in the presence of dredge-up. For exam- 
ple, let us suppose that a superwind suddenly ejects all 
the envelope of the 3M Q star at L = 15000 L Q . Figure [l3] 
then predicts that the mass of the white dwarf's remnant 
would be ~ 0.66 M in the presence of dredge-ups. This 
is much lower than the predicted white dwarf's mass of 
0.73 M Q predicted in the absence of dredge-up. Mass loss 
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Fig. 13. Same as Fig. [L2|, but for the surface C/O ratio 
as a function of surface luminosity at maximum pulse ex- 
tension (solid line) or at the luminosity dip during the 
interpulse phase (taken here to be half of the luminos- 
ity at maximum pulse extension, dashed line). The dotted 
line divides the regions where C stars (C/0>1) and M 
stars (C/0<1) are. The hatched area indicates the range 
of luminosity predicted for given C/O ratios. 

combined with dredge-ups thus further helps to constrain 
the mass of white dwarf remnants. 

7.3. Carbon star formation 

Let us now follow through our synthetic calculations the 
carbon abundance predicted at the surface of a 3 M Q star. 
The amount of carbon dredged- up to the surface, M^ up , 
is given by 



Mg p = X p ("C) x M dup) 



(19) 



where X p ( 12 C) is the 12 C mass fraction left over by the 
pulse in the intershell layers. From the values reported in 
Table |, a mean value of X P ( 12 C) = 0.21 can be assumed. 
The initial 12 C and 16 mass fractions in the envelope are 
taken from our standard 3 M star at the beginning of the 
TP-AGB phase, and are equal to 2 x 10~ 3 and 9.7 x 10~ 3 , 
respectively. The resulting evolution of the surface C/O 
number ratio is shown in Fig. It is seen that a 3 M 
star satisfying Eq. [lo] becomes a C star after about twenty 
dredge-up episodes, at L ~ 15000 L . Furthermore, it is 
known that the luminosity decreases by a factor of almost 
two during about 20% of the interpulse phase. This means 
that the star displayed in Fig. [l3| could be observed as a 
C-star at luminosities as low as 7500 L . 

I should stress that the effects of mass loss have been 
neglected in those synthetic calculations, being outside the 
scope of this article. Mass loss would decrease the dilution 
factor of 12 C into the envelope, but would probably also 
decrease the dredge-up efficiency. Those effects should be 
considered in a future, more detailed, study of C star for- 
mation. 
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8. Conclusions 

The analysis presented in this paper addresses several 
questions regarding the occurrence of 3DUP in AGB star 
models, the laws governing that phenomenon, and some of 
its implications on the structural and chemical evolution 
of those stars. The analysis is based on model calculations 
of a 3 M star using the Schwarzschild criterion and per- 
formed cither with or without extra-mixing. 



Modeling dredge-up in AGB stars The use of a local 
criterion to delimit convection borders, such as that of 
Schwarzschild, and without any extra-mixing is shown to 
lead to an unphysical situation which prevents the deep- 
ening of the convective envelope into the H-depleted re- 
gions, and which thus also prevents the occurrence of the 
3DUP (Sect. |). The Schwarzschild layer is undefined in 
those models due to the development of a discontinuity 
in the hydrogen abundance profile at the bottom of the 
convective envelope. Model calculations using no extra- 
mixing confirm the failure to obtain dredge-up in those 
conditions (Sect. ||). That conclusion is found to be inde- 
pendent of the mixing length parameter, stellar mass, or 
numerical space and time resolution of the models. Mod- 
els using the Schwarzschild criterion with no extra-mixing 
are thus inadequate to describe the 3DUP phenomenon. 

Models of the same 3 M Q star but using extra-mixing, 
on the other hand, lead to efficient dredge-ups from the 
11th pulse on (Sect. ||). This results directly from the un- 
stable character of the lower boundary of the envelope 



against extra-mixing (Sect. 2.2). The calculations further 
reveal that the dredge-up characteristics are insensitive 
to the extra-mixing parameters (such as the extra-mixing 
extent and efficiency). This important conclusion is most 
welcomed since nothing is known yet about the character- 
istics of the extra-mixing expected to occur in AGB stars. 
It enables, in particular, to study the dredge-up character- 
istics with some confidence and without worrying about 
the extra-mixing parameters. 

Although the dredge-up predictions are rather insen- 
sitive to the extra-mixing parameters, a proper handle of 
the extra-mixing procedure is essential in order to obtain 
reproducible predictions. In particular, a purely numeri- 
cal extra-mixing can lead to model-dependent predictions 
(Sect. 5.3; see also Mowlavi 1999a). The use of instan- 
taneous mixing in the extra-mixing region can also lead 
to convergence difficulties. The calculations presented in 
this paper use a diffusive overshooting with a decreasing 
bubble velocity field in the extra-mixing region. It leads 
to smooth chemical abundance profiles and to a proper 
location of the Schwarzschild layer. 

The dredge-up process The insensitivity of the dredge-up 
predictions on the extra-mixing parameters results from 
the fact that the dredge-up rate is limited by the time scale 



of the thermal readjustment of the envelope (Sect. p.l| ). As 
dredge-up proceeds, H-depleted matter is lifted from the 
core into the envelope and its thermal state must relax to 
that of the envelope (see also Mowlavi 1999b). The dredge- 
up rate is estimated to 1CP 5 — 10~ 4 M Q /yr from model 
calculations (Sect. |5|). This is also what is expected from 
an analytical estimation using simplified arguments of the 
physics involved in the dredge-up process (Sect. |6.1| ). 

Dredge-up laws The model calculations performed with 
extra-mixing on selected afterpulses of the standard 3 M© 
star predict linear relations both between the dredge-up 
rate and core mass (Eq. and between dredge-up effi- 
ciency and core mass (Eq. |10| ) . In those calculations where 
the feedback of dredge-up on the AGB evolution is not 
taken into account, the dredge- up rate and efficiency fur- 
ther depend linearly on the stellar luminosity. 

When the feedback of dredge-up on the AGB evolution 
is taken into account, the analysis presented in Sect. || sug- 
gests that the dredge-up rate keeps its linear dependence 
on the core mass. This conclusion is valid if the M c —R c —L 
relation is of the form of Eq. [l4|, R c being the radius of the 
H-depleted core. In that case, we expect the dredge-up ef- 
ficiency to level off at A = 1, at which point the core mass 
docs not increase anymore from one pulse to the next. The 
asymptotic core mass of the 3 M Q star model presented in 



this paper is M* = 0.71 M Q (Sect. |7.l|) 



Consequences on the AGB evolution The synthetic calcu- 
lations performed in Sect. [?] confirm the asymptotic evo- 
lution of the dredge-up efficiency towards unity and of the 
core mass towards M* . This helps to constraint the initial- 
final mass relation of white dwarfs, M* giving an upper 
limit for the mass f the white dwarf remnant. 

The surface luminosity is not as much affected by the 
dredge-ups as is M c . This results from the adopted M c — 
R c — L relation combined with the continuous decrease of 
R c . The non- validity of the M c —L linear relation, assumed 
in Eq. |l8|, is made evident in Fig. [l^. The increase in L 
with time predicted by the synthetic calculations is only 
slightly less steep than the L(t) increase predicted in the 
standard calculations without dredge- up (Fig. [Tl"| ). 

Mass loss can further constrain the initial-final mass 
relation of white dwarfs, because it depends on the lumi- 
nosity through the stellar radius and effective temperature 
(Bowen & Willson 1991), and because the luminosity in- 
creases with time despite the leveling of M c . We thus ex- 
pect mass loss [or superwind(s)] to peal off the AGB's en- 
velope at a lower core mass than in models not experienc- 
ing dredge-up. For example, if a superwind is assumed to 
eject all the envelope when the star reaches L = 15000 L Q , 
the 3 Mq star models with dredge-up predict the forma- 
tion of a white dwarf remnant of 0.66 M Q . This has to be 
compared with a predicted remnant of 0.73 Mq in the ab- 
sence of dredge-up. The synthetic calculations presented 
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in Sect. [7] are, however, preliminary and based on simpli- 
fied assumptions. Future investigations along these lines 
should refine the quantitative predictions given in this pa- 
per. 

Finally, a 3 M Q carbon star is predicted to form after 
about twenty pulses experiencing dredge-up (i.e. at about 
the 25th-30th pulse on the AGB), at L ~ 15000 L . Tak- 
ing into account the luminosity dip experienced during 
about 20% of the following interpulse time, such a C star 
could be observed at luminosities as low as 7500 L . 
Acknowledgments. I thank Drs. A. Jorissen and G. Meynet for 
their careful reading of the first version of this article, and Dr. 
G. Meynet for many useful discussions. 
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